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Abstract 

A method of calculating Feynman diagrams from their small momentum 
expansion |Q] is extended to diagrams with zero mass thresholds. We start 
from the asymptotic expansion in large masses Q (applied to the case when all 
Mf are large compared to all momenta squared). Using dimensional regular- 
ization, a finite result is obtained in terms of powers of logarithms (describing 
the zero-threshold singularity) times power series in the momentum squared. 
Surprisingly, these latter ones represent functions, which not only have the 
expected physical "second threshold" but have a branchcut singularity as well 
below threshold at a mirror position. These can be understood as pseudothresh- 
olds corresponding to solutions of the Landau equations. In the spacelike region 
the imaginary parts from the various contributions cancel. For the two-loop 
examples with one mass M, in the timelike region for ~ we obtain ap- 
proximations of high precision. This will be of relevance in particular for the 
calculation of the decay Z —> hb \n the mb = approximation. 
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1 Introduction 



Once it has been observed [Q] that the calculation of Feynman diagrams on 
their cut can be performed with high precision from their Taylor expansion 
coefficients, there are several advantages of this method, which make it really 
quite attractive: firstly, the Taylor coefficients being known, the remaining cal- 
culation of the diagram in the whole complex plane is a relatively easy task. 
Secondly, more important, the precision with which the coefficients can be 
calculated (from vacuum diagrams) is practically unlimited (e.g. 50 to 100 
decimals with the multiple precision program of is "standard"). This last 
property is of particular relevance in higher loop orders when many diagrams 
(of the order of 1000) contribute, namely the high precision of the Taylor coef- 
ficients suggests that in such a case the scalar amplitudes should be added on 
the level of their Taylor coefficients. Finally, as mentioned in Q, in the 2-loop 
3-point case, in which we are mainly interested here, there occur in general 10 
numerator scalar products of 4 momenta, but only 9 (internal or external) lines 
against which to cancel these. This causes serious problems in the evaluation 
of two-loop vertex Feynman diagrams which are not present if only bubble in- 
tegrals are to be evaluated like in the Taylor expansion. For these reasons it 
appears worthwhile to develop and extend the method of Taylor expansion fur- 
ther to make it applicable to the various kinematical situations in the Standard 
Model. 

The purpose of the present paper is to demonstrate an extension of the 
previous method Q, which can be applied for vertex diagrams with massless 
thresholds. In such a case the Taylor expansion of the Feynman diagrams does 
not exist because of logarithmic singularities at zero momenta squared. The 
method which we propose for this type of diagrams is a combination of using 
standard explicit formulae for asymptotic expansions in large masses P| (see 
[Ql for a short review) and the summation procedure of [Q]. Thus, in the large 
mass limit (M^ ^| |) we get power series in /M"^ factorizing powers of 
ln(-gV^^)- These power series can be summed by means of Fade approximants 
such that the validity of this expansion is extended to large g^/Af^, in the 
spacelike region as well as in the timelike region (due to conformal mapping). 
Note that the general formulae for asymptotic expansions in momenta and 
masses have been successfully applied in a number of papers [^, ^- In 
particular, two-loop self-energy diagrams with general masses were calculated 
in in the region of small and large momentum, and in ||^] in the case of 
massless thresholds and thresholds with small masses. 

Our paper is organized as follows: in Sect. 2 an introductory example of a 
two-loop self-energy diagram is given for which the factorization of the 
ln(-(jf^/M'^) is explicitly known. In Sect. 3 we show cases of interest for the 
decay Z ^ bb and select typical examples for the demonstration of our method. 
Sect. 4 recalls the general method and demonstrates the calculation of the 
"naive" part. In Sects. 5 and 6 the above examples are worked out explicitly 
and Sect. 7 contains our conclusions. 
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2 An example of a two-loop self-energy 
diagram with zero threshold 

Before turning to complicated vertex diagrams we demonstrate how the sum- 
mation by Pade approximants works in a simpler case of the self-energy integral 



Js (see Fig. 1) with zero threshold |10|, where an explicit result written in the 



form of the large mass expansion is known: 




with 




3 ^'^(i) ( __y \ (2) 

1 y n!r(i) , , , 




2M2^r(n + |)(n + l 



Summing the series for /o and /i in the "standard" manner Q, the results of 
Table 1. are obtained. They demonstrate that high precision numerical values 
can be obtained for a large range of g2 values in this way. 



3 Two-loop vertex diagrams with zero 
thresholds 

Concerning the vertex diagrams, there are many different topologies contribut- 
ing to a 3-point function in the Standard Model. For our purpose of demon- 
strating the method, we confine ourselves to the "planar" case with the topology 
shown in Fig. 2 (see "generic"). Distributing in all possible ways massive par- 
ticles over the six virtual lines of this figure, one finds 22 diagrams with a zero 
mass threshold. Our main interest in the present study is, however, besides 
the development of a new method, its possible application to the calculation 
of the decay amplitude for Z bb. For this process mb = can be consid- 
ered as a good approximation and for the corresponding kinematical situation 
{Pi = pi = 0) the new approach is directly applicable (see also [Q]). Therefore, 
instead of presenting all possible 22 diagrams with zero mass thresholds, in 
Figs. 2 and 3 we give only those contributing to Z — > 66 (with mf, = and 
nit large, neglecting quark mixing). Fig. 3 presents infrared divergent diagrams 
and they are merely given for completeness: Cases 7,8 and 9 will be consid- 
ered in a separate publication. Case 10 is a diagram with massless particles 
only, and our method is of no relevance here. This diagram (and also other 
topologies), however, has been evaluated in ||8| (see also |^). In Figs. 2 and 3 
we choose all non-zero masses to be equal (= M). 
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Typical examples which we work out here explicitly are Case 1 and Case 
5 of Fig. 2. These have in the expansion coefficients terms of the order 1/e 
{Case 1) and 1/e and {Case 5), respectively (c? = 4 — 2e). Taking into 
account factors of the form (^jj^^ and il = Pi ~P2): (see below) and 

expanding these in terms of e, the poles in e as well as the dependence on the 
scale parameter (j. drop out. In the final result there remain terms factorizing 
ln{—q^/M'^) and ?n^(— g^/M^) from the latter expansion after cancellation of 
the corresponding powers of e. 

4 Large mass expansion and calculation of 
the "naive" part 

The large mass limit is obtained in the following manner: if some masses are 
much greater than the other masses (in our case all the small masses are zero) 
and all the momenta, one has |Q] 

M — !■ oo . , 

Fr{pi,P2,M;e) ^ Y.PrhiPuP2;e) o%.F^{q\M;e), (4) 

7 

F.. standing for (sub-)diagrams and reduced diagrams characterized by their 
index: F for the original diagram, 7 a subdiagram (7 C F) and F/7 obtained 
from the original diagram by factorizing the product of scalar propagators as 
Ilr = Ilr/^n^ such that more explicitly we have {I is the number of loops) 

F^/^f oTqiF^ = J dki ■ ■ ■ dkiU-p/^TgilL^. (5) 

Here T . is the operator of the Taylor expansion w. r. t. the set of external 
momenta q'^ of the subgraph 7. The summation in (Q) is performed over the 
following subgraphs: 

— each 7 (it may be disconnected) contains all the lines with large masses, 

— each 7 is IPI w. r. t. lines with small masses. 

The ("naive") contribution from the original diagram F itself is obviously 
always contained in the J2-y various terms of this sum the integrations of 
the type (^) yield in general divergent coefficients of the asymptotic expan- 
sion. These divergences are both of infrared and ultraviolet nature, the latter 
being due to high powers of integration momenta produced by the T . oper- 
ator. Therefore @) is to be understood in terms of some regularization for 
which we take dimensional regularization. This is so even if the original dia- 
gram is convergent: summing all contributions {^^) the divergent terms and 
those depending on the scale parameter /j, (from dimensional regularization) 
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must cancel. This will be used below as a strong check of our calculational 
procedure. 

In the following we will calculate the diagrams Case 1 and Case 5 of Fig. 2. 
These are typical in the sense that they have one and two zero mass thresholds, 
resulting in the above mentioned terms up to 1/e and in the Taylor coef- 
ficients, respectively or ln{—q^/M'^) and l'n?{—q'^/M'^) in the final result. The 
three particle zero mass threshold in Case 6 does not induce higher than 
terms in the Taylor coefficients either. 

The "higher" terms in the (i.e. all except for the naive one) can be 
handled in a straightforward though tedious calculation. The reason for the 
relative simplicity is that only factorizing massive one-loop bubble (vacuum) 
integrals and massless propagator type integrals occur. Results for these higher 
terms will be given in the following sections. 

The situation is different, however, for the "naive" contribution. In this 
case the approach of [|ll| turns out to be particularly adequate. The general 
expansion of (any loop) scalar 3-point function with its momentum space rep- 
resentation C{pi,p2) can be written as 

oo 

c{pi,P2)= E (^in^nipiyipiripmr, (6) 

l,m,,n=0 

where the coefficients aimn are to be determined from the given diagram. As we 
will show, in the cases under consideration the representation of the Taylor co- 
efficients given in yields only "genuine two-loop bubbles" but no factorizing 
one-loop ones. Introducing the abbreviations (see also Fig. 2) ci = k\—m\,C2 = 
ki — 7712, C3 = ^2 — 771^, C4 = k2 — 7n\ and C5 = /c^ ~ "^ii ce = (^i — k2)'^ — mg, 
we have for the expansion coefficients 



on j- 1 

(i7r2)2aoon = —-r{l^^?' / d''k^d''k2Fn , (7) 

71 + 1 J Ci C2 C3 C4 C5 C6 



on 1 

Ci C2 C3 C4 C5 C6 

where is the scale parameter of dimensional regularization. 

In general the vertex diagram depends on three external momenta squared 
(see ^ , each of which is an independent expansion variable. Putting pf = P2 = 
the corresponding summation indices are also zero. In (j^ Fn is given by 



2 ~ "T-eJ 



' i^i + ^1 - mlr+'''-^^'-'^{2kik2r-^ • C6 J , (8) 

0=1, odd 



the coefficients a"^, being known explicitly for arbitrary d [11 1. We see that 
due to cancellation of cq in the above sum over a this contribution contains 
only factorizing one-loop terms which vanish for the above mass combinations 
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in dimensional regularization. Thus it remains to calculate the "genuine two- 
loop" contributions. 

In Case 1 (me = M) we write {X = u + u' — 2ji') 



13=0 y 



and obtain 



aoon 



.2\2e 



x{MY~^^^^ I Ad'^k2jr^ (9) 

with = n - (i/ + v') + /j/ + /3 + 2 and V2 = {i^ + v') - (3 - ^ + 3. The two-loop 
bubble integrals can be evaluated explicitly: 



^ f.,+.,+i) r(^i + + 1 - t^)r(i - + i)r(^i + 1 - j) 
^ ' r(z/2)r(f)r(2i/i + i/2 + i-d)(M2)('^i+^2+i-<i) 

Expanding in £, a divergent ^ - term is obtained, which in this case comes from 
the infrared divergence of the above integral. 

Similarly as above we calculate the "naive" part for Case 5: with 

{kir'iki + ki - M-Y = E E (^) (^) {-My-^ikif-^ikiY'^-^ 

we have 



2" ^ V 1 4^4^(>\(f>\ 
><(-^V-|$/d^*.<i^fe(^lpr^j^ (10) 

with ui and U2 as above. The situation is now somewhat more complicated 
due to the fact that 777.3 / f^b and the following partial fraction decomposition 
needs to be performed ( p = i/2 — 2, 7773 = and 7775 = M): 

1 ^ 1 1 1 1 

■ ^ /'A#2\«+l-i J+1 + 
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yielding the two-loop bubble Integrals 



If 1 

( + i + 2 - d)T{l - i - i)r(^i + i + 1 - f )r(f - u,) 
^ ' r(z^i)r(f)r(i + i)(M2)K+i+2-'^) 

and 



^ ^ r(f)r(2z^i + 2-(i)(M2)('^i+2-d) 

Expanding in e, in this case ^-terms occur. 



5 Threshold singularity of the type 

ln(-gVM2) (Case 1) 

According to Fig. 2 there is only one subdiagram 7 which contains all four 
massive lines, namely the box. Accordingly the only contribution with 7 7^ F 
reads 



1 

[(fci - ^2)2 - M2][(g + A:2)2 - M2][(pi - A:2)2 - M2] ■ ^^^^ 

Here the external momenta for the box subgraph are besides pi and q also 
the loop momentum ki. Note that, for the contributions of each subgraph 7 
in one can choose loop momenta in an appropriate way. For example, in 
( [TTI) the rooting of the loop momenta is different from the "naive" contribution 
(observe also below similar rerooting in Case 5). 

The evaluation of (11) results in products of one- loop massive bubble inte- 



grals with monomials in the numerator and one-loop propagator-type massless 
integrals, the latter yielding the factor (—1/q'^Y. Explicitly we obtain 

' ' %\i-i\f:c'^\^)i'i'/Mr, (12) 



{4ttY M4 I M2 



n=0 



where 



.mU)= V \ i)(h+i2+m)/2 jn - ih +i2- n3)/2)! 

H.2.n3>0.lTt.2+n3 even .fe) (n - (zi + ^2 + ns)/2)l 

*1 +^2 "'""'3 
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iili2^.0{ii + 12 - - k + jsWj-h + ^2 + js) 

{{ns - J3)/2)!((ii +12- i3)/2)!((ii -12+ j3)/2)!((-ii +i2 + j3)/2)! 

^^ r(e)r(l-e)r((ii+Z2-j3)/2 + l-6) 
r((ii+i2-i3)/2 + 2-2e) 
xC(n + (ii + i2 + n3)/2 + 4; (ii + is + n^)/!) (13) 

and 

C(r,s)=r(^r-s-0/r(r). (14) 

The e— poles in this case come from ultraviolet divergences. Adding the 
naive and the above contribution yields (up to a factor (l/lGvr^)^) 

( 00 00 
F^(g^M2) = ^ E foniqVMT + E /^(gVM^)- ln(-gVM2) 

l?i=0 n=0 

= ^ {/o(gVM2) + fi{q'/M') H-qVm')} , (15) 

i.e. all e— poles and the scale parameter have cancelled, which at the same 
time serves as a helpful check of the correctness of the calculation. Evaluating 
the complicated coefficients and summing all contributions was performed with 
FORM [|l^]. For the expansion coefficients of /i we found by inspection the 
following recurrence relation: 



8(2n + 5)(2n + 7)(2n + 9)(n + 4) /i,„+3 = 
-4(2n + 5)(2n + 7)(n + 4)2 /i^„+2 

+ 2(2n + 5)(n + 3)(n + 2)2 /i,n+i + (n + 3)(n + 2)3 , (16) 
from which we obtain: 

r(n + 2)r(|) ^ i + 2(-i)^- 

6(-4)«r(n+|)^t^, {j + iy ■ ^''^ 

This explicit form allows to sum the series for /i, yielding the following 
integral representation: 



n=0 



2;, 



M2 



(18) 



where 



(z,x) 



4 ar csin ( y/zx/i) 



(1 — z)y/zx{4: + zx) 



(19) 
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An interesting feature of this representation is that two thresholds, one at 
= 4M^ and another at = — 4M^ can immediately be read off. There 
ought to be, however, only one threshold at = 4M^. Indeed it turns out 
that /o has a singularity at = — 4M^ as well and that for q^ < — 4M^ the 
imaginary part of /i ln(— g^/M^) and /o cancel. 

Note that the position of the "mirror" threshold of the functions /o and fi 
exactly corresponds to a pseudothreshold of the given Feynman diagram which 
is a solution of the Landau equations ( see e.g. [0]). In fact solving these, we 
obtain for rrii = m2 = the pseudothreshold = — (ms + mg)^. ^ An even 
simpler case is the corresponding one- loop vertex with a massive line connecting 
two massless ones. This graph is proportional to ( x=(7^/M^ ) 



Li2{l + x) - C(2) = -Li2{-x) - ln(l + x) ln{-x), (20) 

where in the second form of writing the logarithmic singularity at x=0 has 
been isolated and now again the pseudothreshold singularity at x=-l appears 
in the separate terms. For two- loop self-energy diagrams one obtains by in- 
spection corresponding results from Refs. 0] and |14|. At this point it is 
interesting to note that the self-energy diagram I^, dealt with in Sect. 2, has no 
pseudothreshold, i.e. no such solution of the Landau equations exists. 

We have found, however, a closed form as in ( [l^ ) only for fi and not for 
ft). Therefore, to demonstrate the above statement, we have to rely on our 
numerical approach. Results for Case 1 are given in Tables 2 to 4 (in all tables 
the results are given up to a factor (jfoF^plvF)' 

Table 2 gives results for spacelike q"^. For > q^ > — 4M^ we achieve fast 
convergence, calculating the indicated Fade approximants, taking into account 
Taylor coefficients up to n=30. Of course the integral is real in this domain and 
agrees excellently with the Monte Carlo control calculation (five dimensional 
integration over Feynman parameters with an estimated relative error < 10~^). 

For < — 4M^ we first performed a conformal mapping in terms of an "w- 
transform" ||l| before using Fade approximants to sum the series for /o and /i. 
This is of course always necessary if one wants to calculate a function on a cut. 
As a result we obtain in the above domain small imaginary parts, which are 
due to insufficient cancellation of the two imaginary parts. It is seen, however, 
that these become quickly small with increasing order of the approximants. 
The situation becomes even more amazing if one looks at the imaginary parts 
of /o and /i separately: they do increase with q'^ relative to their real parts as 
is demonstrated in Table 4! Nevertheless, in this manner we obtain even for 
large negative values like g^/M^ = — 50 an accuracy of at least 3 decimals as is 
also confirmed by the MC control calculation. 

Table 3 gives our results in the timelike region. The analytic continuation 
of the logarithm requires in this case to write 

ln{-q^/M'^) = In \ q^/M'^ \ -m (21) 



^We are grateful to J.B. Tausk who has drawn our attention to this property and provided several 
examples, including the following one. 
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so that an imaginary part is obtained for > 0. Table 3 gives only results 
above the second threshold since below the convergence is even better: at 
= we obtain a precision of 9 decimals with merely 9 coefficients ([4/4]). 
A precision of at least 3 decimals is achieved both for the real and imaginary 
part up to gV^^=50. This can be concluded from the convergence properties 
of the approximants. 



6 Threshold singularity of the type 

ln^{-q^lM'^) {Case 5) 

There are two massless cuts so that we shall have the double logarithm in the 
expansion. The set of subgraphs in this case is given by 71 = F and the higher 
terms with 7 7^ T : 72 = {3456}, 73 = {1256}, and 74 = {56}. Note that 73 
and 74 are disconnected. The subgraph 71 was discussed in Sect. 4. In terms 
of integrals of the form dH) we have 
for 72 



2 



d"fc2^2^p„,,fc, ^(^^ ^ _ _ ^^)2 _ M^){q - ^2)2' ^^^^ 



for 73 



1 ^ 1 

'A;2(g-A;2)2 ^'•'^'(pi + fe2)2-M2 
^ / _i.\2_ ^/riM r, _ h-.\2 (2"^) 



and finally for 74 



The results for these three integrals look as follows. The subgraph 72 yields 

-, -, / \ ^ / o \ ^ rsT, 

1 I jl'^ \ IJL 



Th] E4^^W('^VmY, (25) 



(47r)'iM4 \M^^ ^ J n=0 
where 

c(2)fe)= V i\{ii+i2+m)n - (^1 + ^2 - ^3)72)! 

n.2."3>0.lTt«2+n3 even (n - (zi + ^2 + n3)/2)! 

+12 +ri2 <2n 
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iili2^.0{ii + 12 - - k + jsWj-h + ^2 + js) 

{{ns - j3)/2)!((ii +12- i3)/2)!((ii -12+ j3)/2)!((-ii +i2 + j3)/2)! 

^^ r(e)r(l-e)r((ii+Z2-j3)/2 + l-6) 
r((ii+i2-i3)/2 + 2-2e) 
xC(2 + ii + is, 2 + n - (ii + ^2 - n3)/2; (ii +12 + n3)/2),(26) 



and 

C(n,r2; s) = ' '^.J,^;^^:^^, ' ^^ (27) 



r(ri + r 2-g-f)r(5-r2 + ^. 

r(ri)r(.+ 2. 



Note that this expression is obtained from the corresponding contribution of 
the non-naive part (O) for our first diagram by the change 



C(2 + ii+i2,2 + n- {ii+i2- nz)/2- {h + i2 + nz)/2) 
C{n + (n +i2+ n^)/2 + 4; (n +12 + n3)/2). 

The contribution from 73 takes the foUowing exphcit form: 
where 

ii ("1,-2) 

c^f^e) = (-!)" E EE (-1) 



n+ii. 



H-("l-Jl)/2<" 

r(2 - f )r(f - i)r(n - n + (m - ji)/2 + f - 1) 

r(n-ii + (ni-ji)/2 + f) 

xC(ii + l,ii + 2;(ni+ji)/2). (29) 



Finahy, for 74 we obtain 

1 1 ( fJ- 



2X 2e 



(J^5^(-^j 54''(^)hVM'r, (30) 
with 

After summing up ah four contributions we see that the double and single 
poles in e cancel as well as the scale parameter with the result 

I 00 2 

Fr{q',M') = E E/."(«VM^rin^(-'zVM2) 
= ^ {foiqVM^) + h{q^/M^) H-q^/W) + ^(gVM^) \n\-q^/M'')] , (32) 
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where the fjn are now obtained in terms of rational numbers and a ("(2) con- 
tained in /on- 

The next surprise is that f2{x) can be summed analytically, yielding 
(x = q^M^) 

72(2;) = /n^(l + x)/x^ for X > —1, 

= [In I 1 + X I — ivr f/x^ for X < -1. (33) 

Thus, as in Case 1, we have to Pade approximate /o and /i only. The C(2) 
term in /o can in principle also be summed, i.e. 

/o(x) = /o(x) + 2C(2)/2(x), (34) 

with fo{x) having only rational numbers as expansion coefficients. This split- 
ting does, however, worsen the convergence properties of the series for /o and 
therefore has no practical meaning. Apart from that the situation is similar as 
in Case 1, i.e. /o has two thresholds starting at = ztM^, respectively (in- 
deed solving the Landau equations with mi = m2 = ms = m4 = one obtains 
the additional pseudothreshold = —rri^ ). The function /i, however, is real 
for timelike q^. In the spacelike region the cancellation between the imaginary 
parts has now to take place between three complex functions. This is shown in 
Table 5, but apparently this cancellation is not as perfect as in Case 1. This is 
due to a relatively bad convergence of the Pade's for /o (with 30 coefficients) 
while the convergence properties of /i are much better and /2 is anyway given 
analytically. Fortunately the convergence of /o is much better in the timelike 
region (see Table 6) so that up to q'^ = lOM^ the achieved precision is at least 3 
decimals. Close to the second threshold at q'^ = the convergence is indeed 
excellent. It should be noted that for the physical application we have in mind, 
i.e. Z — >■ bb, this is just the case of interest. It is worthwhile to note the sharp 
increase for low q"^, in particular in the timelike region, due to ln'^{-q^/M^). 

7 Conclusions 

We have presented a modification of the method of Q for the case when massless 
thresholds are involved. To do this, the starting point was the asymptotic large 
mass expansion, rather than the Taylor expansion in the external momenta. 
Using the explicit formulae for the coefficients of this expansion we obtained a 
finite sum of powers of the logarithm of the external momenta times power series 
in q^ /M"^ and then applied the technique of conformal mapping and summation 
by Pade approximants to each of these series separately. 

We have shown that this new strategy enables us to obtain high precision 
numerical values in the domain of physical interest including the region beyond 
the second threshold in spite of the fact that the initial Feynman diagrams 
have their first threshold at the origin of the complex plane. Below the second 
threshold, in a certain domain around q^ = 0, the functions factorizing powers 
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of ln[—q^/M'^) can in general be assumed to be analytic so that even without 
Pade approximants their power series converge. At the second threshold with 
definiteness one can only say that at least one of them must have a cut starting. 
This is confirmed by our calculation, i.e. the application of conformal mapping 
and Pade approximants (recall the observation that in Case 5 /i is real for 
timelike q^). The occurrence of "ghost thresholds" appears to be related to 
pseudothresholds obtained from the Landau equations. 
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Figure 1: Self-energy integral I3; dashed lines are massless and solid lines massive. 
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Figure 2: Planar diagrams with zero thresholds occurring in Z — ^ hh. 
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Figure 3: On-shell infrared divergent planar diagrams with zero thresholds. 



Table 1. Results for the propagator diagram I3 
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Table 2. Results for spacelike (Case 1) 
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Table 3. Results for timelike (Case 1) 
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Table 4. Functions /o and /i for spacelike q^ (Case 1) 
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Table 5. Results for spacelike (Case 5) 
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Table 6. Results for timelike (Case 5) 
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